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l/-) '. ABSTRACT 

I/"") ' Galaxy-galaxy lensing is rapidly becoming one of the most promising means to accu- 

' rately measure the average relation between galaxy properties and halo mass. In order 

to obtain a signal of sufficient signal-to-noise, one needs to stack many lens galaxies 
according to their property of interest, such as luminosity or stellar mass. Since such 
a stack consists of both central and satellite galaxies, which contribute very different 
f~| . lensing signals, the resulting shear measurements can be difficult to interpret. In the 

Q^' past, galaxy-galaxy lensing studies have either completely ignored this problem, have 

applied rough isolation criteria in an attempt to preferentially select 'central' galax- 
ies, or have tried to model the contribution of satellites explicitely. However, if one is 
^ ' able to a priori split the galaxy population in central and satellite galaxies, one can 

measure their lensing signals separately. This not only allows a much cleaner measure- 
ment of the relation between halo mass and their galaxy populations, but also allows 
a direct measurement of the sub-halo masses around satellite galaxies. In this paper, 
we use a realistic mock galaxy redshift survey to show that galaxy groups, properly 
' selected from large galaxy surveys, can be used to accurately split the galaxy popula- 

tion in centrals and satellites. Stacking the resulting centrals according to their group 
mass, estimated from the total group luminosity, allows a remarkably accurate recov- 
ery of the masses and density profiles of their host haloes. In addition, stacking the 
corresponding satellite galaxies according to their projected distance from the group 
center yields a lensing signal that can be used to accurate measure the masses of both 
sub-haloes and host haloes. We conclude that an application of galaxy-galaxy lensing 
measurements to group catalogues extracted from large galaxy redshift surveys offers 
a unique opportunity to accurately constrain the galaxy-dark matter connection. 
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1 INTRODUCTION 

Understanding the connection between galaxies and dark 
matter haloes is a major challenge in modern astrophysics. 
From the perspective of cosmology, the galaxy-dark mat- 
ter connection is required in order to translate observations 
of galaxy clustering in terms of the distribution of (dark) 
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matter. From the perspective of galaxy formation, it reveals 
how halo mass impacts on the properties of the galaxies 
that form within them, and thus how the various phys- 
ical processes that play a role in galaxy formation scale 
with halo mass. Consequently, a great amount of effort has 
been devoted to establishing the galaxy-halo connection, ei- 
ther through numerical simulations (e.g., Katz, Weinberg & 
Hernquist 1996; Fardal et al. 2001; Kay et al. 2002; Springel 
et al. 2005; Springel 2005), through semi-analytical model- 
ing (e.g., White & Frenk 1991; Kauffmann et al. 1993, 2004; 
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Mo & Fukugita 1996; Mo, Mao & White 1999; Somerville & 
Primack 1999; Cole et al. 2000; Benson et al. 2002; van den 
Bosch 2002; Croton et al. 2006), or via statistical approaches 
(e.g., Jing et al. 1998; Seljak 2000; White 2001; Berlind & 
Weinberg 2002; Yang et al. 2003b; van den Bosch et al. 2003; 
2005c; Zheng et al. 2005; Cooray 2005, 2006). However, nei- 
ther of these methods provides a direct diagnostic of the 
galaxy-halo connection. 

Direct measures of the dark matter haloes around galax- 
ies come either from dynamical tracers or from gravitational 
lensing. Galaxy rotation curves and strong lensing, although 
extremely powerful, typically only probe the inner part of 
the dark matter halo, and are therefore not well suited to 
determine the total halo mass. The only dynamical trac- 
ers that probe the gravitational potential sufficiently far out 
to allow for an accurate mass estimate are satellite galax- 
ies (Zaritsky & White 1994; McKay et al. 2002; Prada et 
al. 2003; Brainerd & Specian 2003; van den Bosch et al. 2004; 
Conroy et al. 2005). However, arguably the best method 
to directly probe the dark matter haloes around galaxies 
is galaxy-galaxy lensing, whereby the gravitational field of 
lensing galaxies induces small tangential shear distortions in 
the images of distant background galaxies (e.g. Natarajan 
& Kneib 1997). Unfortunately, since the weak lensing sig- 
nal around individual galaxies is too small to be detected, 
one can only infer ensemble-averaged properties: by combin- 
ing the signal from a large number of lensing galaxies one 
obtains the galaxy-mass cross correlation function, which 
can be used to infer the mass distribution around galaxies, 
but only in a statistical sense. Note, however, that a similar 
problem hampers the satellite-dynamics method, where the 
number of satellite galaxies of individual host galaxies is too 
small to allow for a reliable estimate of the halo mass, and 
ensemble averaging has to be used as well. 

Since the galaxy-galaxy lensing signal is so weak, and 
one has to carefully correct for a number of observational 
effects, such as anisotropies in the point-spread function and 
shear induced by the camera optics, it took twelve years 
since the first attempt by Tyson et al. (1984) until the first 
detections of galaxy-galaxy lensing (Brainerd et al. 1996; 
Griffiths et al. 1996; dell'Antonio & Tyson 1996). In recent 
years, however, the progress has been enormous, largely due 
to the advent of deep galaxy surveys with large sky coverage 
(e.g. Fischer et al. 2000; Wilson et al. 2001; Smith et al. 2001; 
Mckay et al. 2001; Hoekstra et al. 2003, 2004, 2005; Sheldon 
et al. 2004; Kleinheinrich et al. 2006; Parker et al. 2005; 
Mandelbaum et al. 2006a,b). 

The fact that galaxy-galaxy lensing only yields en- 
semble averaged properties complicates the interpretation. 
Since galaxies with different luminosities and with different 
morphologies are expected to reside in haloes of different 
masses, the ensemble averages are complicated sums over 
large ranges in halo mass. Unless one has a prior knowledge 
of the relation between galaxy properties and halo mass, an 
unbiased interpretation of the measurements is basically im- 
possible. For example, in the past most studies have used the 
observed lensing signal only to constrain the normalization 
of an assumed relation between galaxy luminosity, L, and 
halo mass, M (e.g., Fisher et al. 2000; Smith et al. 2001; 
Wilson et al. 2001; Hoekstra et al. 2003). When redshift in- 
formation regarding the lenses and/or sources is available, 
one can somewhat improve the constraints. For example, 



Hudson et al. (1998) were able to simultaneously constrain 
the normalization and the slope of the L-M relation, due 
to the fact that they had photometric redshifts available 
for both their lenses and their sources. Nevertheless, their 
results clearly depend on the assumed functional form (a 
power-law) of the L-M relation. 

Therefore, one of the ultimate goals in galaxy-galaxy 
lensing is to be able to measure the lensing signals for galax- 
ies with different intrinsic properties. This not only requires 
a very large sample of galaxies (both lenses and sources), but 
also redshifts for the lensing galaxies, in order to be able to 
(i) convert angular sizes into physical sizes, (ii) compute the 
luminosities of the lenses, and (iii) reduce uncertainties in 
the geometry of the lens-source system. The Sloan Digital 
Sky Survey (SDSS) has proven to be ideally suited to make 
progress along such a direction, and several measurements 
of the galaxy-galaxy lensing signal as function of galaxy lu- 
minosity and morphological type have already been made 
(e.g., McKay et al. 2001; Sheldon et al. 2004; Mandelbaum 
et al. 2006a). 

There is one additional problem that complicates the 
interpretation of these lensing signals, namely the distinc- 
tion between 'central' galaxies, which reside at the center 
of a dark matter halo, and 'satellite' galaxies, which are lo- 
cated on an orbit around a central galaxy. In the current 
paradigm, these satellite galaxies are thought to be asso- 
ciated with dark matter sub-haloes, which are haloes that 
reside and orbit within a larger virialized dark matter halo. 
The idea that satellite galaxies are related to the population 
of dark matter sub-haloes is consistent with both observa- 
tions and numerical simulations (e.g. Kravtsov et al. 2004; 
Natarajan & Springel 2004; Natarajan ct al. 2006; Vale & 
Ostriker 2004, 2006; Conroy et al. 2006; Kang et al. 2005). 
The lensing signal around a satellite galaxy does not only 
reflect the mass distribution of the dark matter sub-halo in 
which it is located, but also that of the larger halo that 
hosts the sub-halo (hereafter host halo). Thus one expects 
the lensing signal from a satellite galaxy to be very different 
from that of a central galaxy (e.g., Appendix B in Yang et 
al. 2003a). In addition, the lensing signal from a satellite 
galaxy not only depends on the masses of the sub-halo and 
host halo, but also on the geometrical orientation of these 
two entities with respect to each other. 

Since a galaxy of a given luminosity and/or morpholog- 
ical type can be either a central galaxy or a satellite galaxy, 
even the lensing signal from a sample of lensing galaxies with 
a narrow range in luminosities and morphological types can 
be difficult to interpret. In the past, galaxy-galaxy lensing 
studies have either completely ignored this problem (e.g., 
McKay et al. 2001; Hoekstra et al. 2003, 2004), have ap- 
plied rough isolation criteria in an attempt to preferentially 
select 'central' galaxies (e.g., Hoekstra et al. 2005; Mandel- 
baum et al. 2006b), or have modeled the contribution of 
satellites explicitely, using either semi-analytical models for 
galaxy formation (Guzik & Seljak 2001; Yang et al. 2003a) 
or a model for the halo occupation statistics (Guzik & Sel- 
jak 2002; Mandelbaum et al. 2005b, 2006a). None of these 
approaches have taken into account the weak lensing signals 
around satellite galaxies at different halo-centric distances. 

In this paper we demonstrate that it is actually pos- 
sible to measure the lensing signal around satellite galaxies 
at different projected halo-centric distances, provided that a 
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well-defined sample of galaxy groups and clusters is available 
to represent galaxy distribution in dark matter haloes. This 
not only allows a much cleaner measurement of the relation 
between halo mass and their galaxy populations, but also 
allows a direct measurement of the sub-halo masses around 
satellite galaxies. Using realistic mock galaxy catalogues we 
show that the halo-based group finder, recently developed 
by Yang et al. (2005a), allows an accurate identification of 
central and satellite galaxies, and that their corresponding 
lensing signals allow an accurate recovery of their mean host 
and sub-halo masses, respectively. Note also that in order to 
perform such a measurement, a large spectroscopic sample is 
needed to group lensing galaxies according to their common 
haloes. Currently, the largest sample of this kind is from 
the Sloan Digital Sky Survey (SDSS). With about 1/3 of 
the eventual SDSS, significant galaxy-galaxy lensing signals 
have been detected separately for three (seven) subsamples 
in galaxy luminosity and colour (e.g. Sheldon et al. 2004; 
Mandelbaum et al. 2006a). This suggests that the present 
SDSS data may already be able to probe the lensing signals 
around galaxies in broad bins of halo mass and halo-centric 
distances. With the completion of the SDSS, and with fu- 
ture deeper imaging surveys in the same part of the sky 
as the SDSS, we anticipate that the effects we are consid- 
ering here can be studied. In a forthcoming paper, we will 
construct shear maps from realistic mock catalogues taking 
into account observational effects to test the feasibility of 
the method we are proposing here. 

This paper is organized as follows. In Section 2 we give 
a detailed description of the weak lensing signal produced by 
host haloes and sub-haloes, separately, in Section 3 we use 
realistic mock catalogues constructed from N-body simula- 
tions to demonstrate how to split the galaxy population in 
central and satellite galaxies, and how their lensing signals 
allow a measurement of their corresponding host and sub- 
halo masses, respectively. Our conclusions are summarized 
in Section 4. 



2 LENSING BY HALOES AND SUB-HALOES 

Galaxy-galaxy lensing measures the profiles of the tangen- 
tial shear, 7t(-R), azimuthally averaged over a thin annulus 
of projected radius R around a set of lens galaxies. This ob- 
servable quantity is related to the mean projected surface 
mass density within the aperture radius R according to 

7«(fl)Ecrit = £(< R) - = AE(JZ) . (1) 

Here E(< R) is the mean surface density within R, E(i?) is 
the azimuthally averaged surface density at R, and 

is the critical density, which is determined by the geometry 
of the lens-source system (Miralda-Escude 1991; see Schnei- 
der 2005 for a detailed review). In the above equation, D\, 
D B and As are the angular diameter distances to the lens, to 
the source and between the lens and the source, respectively. 
Since a uniform mass sheet, such as the mean density of the 
universe, does not contribute to AE, it basically measures 
the excess surface density (hereafter ESD). 



The mean excess surface density around a galaxy is 
specified by the line-of-sight projection of the galaxy-matter 
cross correlation function, £ gm (r), so that 

and 

E(<fl) = -g [ R ydy r ?gm (r)-^L= (4) 
H Jo Jy \/r 2 - y 2 

with p the average background density of the Universe. Note 
in both equations, we have omitted the contribution from 
the mean density of the universe, as it does not contribute to 
the ESD. As we will see below, it is important to distinguish 
between the lensing signal due to host haloes (those haloes 
that are not embedded in a larger virialized structure) and 
sub-haloes (haloes embedded in a host halo). In what follows 
we refer to galaxies at the centers of host- and sub-haloes as 
central galaxies and satellite galaxies, respectively. 

In the halo model, the dark matter distribution con- 
sists entirely of dark matter haloes, and the galaxy-mass 
cross correlation function consists of four terms: three one- 
halo terms, and one two-halo term. The first one-halo term 
is due to the host haloes around central galaxies. Satellite 
galaxies contribute two one-halo terms: one describing the 
density distribution of the dark matter sub-halo and one de- 
scribing the density distribution of the host halo in which 
the sub-halo is embedded. This latter term also depends on 
the relative location of the sub-halo with respect to the cen- 
ter of the parent halo. Finally, the two-halo term describes 
the correlation between the lens galaxy and the large scale 
distribution of dark matter haloes. In this paper we focus 
only on the expected lensing signal at small R, and we will 
therefore ignore this two-halo term in what follows. 

From the above it is clear that central and satel- 
lite galaxies yield different lensing signals (see Hudson et 
al. 1998; Yang et al. 2003a; Guzik & Seljak 2002). In the 
past, various studies have used halo occupation statistics to 
model the total galaxy-mass cross correlation function, due 
to central and satellite galaxies combined (Guzik & Seljak 
2002; Mandelbaum et al. 2005b, 2006a; Yoo et al. 2005). 
Here we adopt a different approach and we investigate the 
AE(i?) of central and satellite galaxies separately. In Sec- 
tion 3 we show that with a decent galaxy group finder one 
can identify central and satellite galaxies, and probe their 
lensing signals separately. 

2.1 Density Distribution of Haloes and 
Sub-Haloes 

As outlined above, the tangential shear due to galaxy-galaxy 
lensing can be used to infer the galaxy-mass cross correlation 
function £ gm (r). On small scales, £g m (r) = p(r)/p — 1, with 
p(r) the azimuthally averaged density distribution around 
the lensing galaxies. This in turn reflects the density distri- 
bution of the host- and of the sub-haloes surrounding central 
and satellite galaxies, respectively. For simplicity, we will 
not treat the baryonic masses of the lens galaxies as sepa- 
rate mass components, but instead consider them included 
in the halo components. As shown in Lin et al. (2006), the 
baryonic component can change the concentration of the to- 
tal mass profile on small scales. We do not expect such effect 
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to have an important impact on our results, because we are 
focusing on lensing signals on scales 7?>50 h~ 1 kpc. 

For the host haloes, we adopt the following density pro- 
file: 



p(r) = 



Po 



(r/r c ) a (l + r/r c ) 3 - 



where 



Po = 



I(c,a) 



dx 



(5) 



(6) 



3/(c,a)' c 3 J Q x a - 2 {l+x)' A - a 

(see e.g., Zhao 1996; Jing & Suto 2000). Note that a = 1 
corresponds to the NFW profile (Navarro, Frenk & White 
1997); a = 1.5 corresponds to the profile proposed by Moore 
ct al. (1999), and a — corresponds to a density distribution 
with a constant density core. Although we will mainly focus 
on haloes with a — 1, we also briefly discuss the impact of 
changing the central cusp slope. This is motivated by the fact 
that (i) observationally the exact value of the inner slope a is 
still uncertain (e.g., Swaters et al. 2003; Simon et al. 2005; 
Sand et al. 2004; Bartelmann & Meneghetti 2004; Gentile 
et al. 2004, 2005; Meneghetti et al. 2005), (ii) numerical 
simulations suggest that haloes may reveal a fair amount of 
scatter in their central cusp slopes (e.g., Jing & Suto 2000; 
Power et al. 2003; Dahle et al. 2003; Navarro et al. 2004), 
and (iii) the density distribution is interpreted as including 
the contribution from the baryons, so that it doesn't have 
to be in perfect accord with that of dark matter haloes. 

For a given value of a, the density distribution of a dark 
matter host halo is specified by two parameters, a char- 
acteristic density po and a characteristic radius r c . Alter- 
natively, one can parameterize the halo by its mass M — 
(4-7r/3)A v j r pr^ ir and concentration parameter c = r v i r /r c . 
Here r vir is the virial radius, defined so that the average 
density within it is A v i r p. Throughout this paper we adopt 
Avi r = 180, and we use Mh and M s to refer to the masses 
of host- and sub-haloes, respectively. 

Numerical simulations have shown that halo concentra- 
tion is correlated with halo mass. Throughout we adopt the 
model of Eke, Navarro & Steinmetz (2001) to model the re- 
lation between c and M, and we assume that it is free of 
scatter. Therefore, the entire density distribution of a dark 
matter (host) halo is completely specified by its mass alone. 

For sub-haloes, we follow Hayashi et al. (2003), who, 
using numerical simulations, found that the density profiles 
of stripped dark matter sub-haloes can be written as 

ft 



Ps{r) = 



■ P(r) ■ 



(7) 



1 + (r/r t ,eH ) 3 

Here ft is a dimensionless measure for the reduction in cen- 
tral density, and r t , c tt is an 'effective' tidal radius that de- 
scribes the outer cutoff imposed by tides. For ft — 1 and 
rt, e s 3> r c , eq. (7) reduces to the original mass profile p(r) 
given in eq. (5), i.e., the mass profile before the sub-halo 
was accreted by the host halo. The parameters ft and rt, e ff 
are determined by the mass fraction of the sub-halo that 
remains bound, f m . Fitting the mass profiles of numerous 
stripped dark matter sub-haloes, Hayashi et al. (2003) ob- 
tained the following relations between r t , e s (in units of the 
characteristic radius of the sub-halo) and f m : 

log[r t ,eff/r c ] = 1.02 + 1.38 log f m + 0.37 (log f m f , (8) 

and between ft and f m : 



log ft = -0.007+0.35 log / m +0.39(log/ m ) 2 +0.23(log/ m ) 3 . (9) 

Using a large cosmological numerical simulation, Gao 
et al. (2004) studied the radial dependence of the retained 
mass fraction f m of a large sample of sub-haloes. Using the 
results shown in their Fig. 15, we obtain the following mean 
relation: 



f m = 0.65(r s /r vir , h ) 2/3 . 



(10) 



where r s is the distance of the sub-halo from the center of the 
host halo, and r v i ri h is the virial radius of the host halo. The 
combination of Eqs. (7) - (10) give a model for the density 
profile of a sub-halo with a given original mass located at a 
given distance from its host halo. 

Note that the above model for the density distribution 
of the host and sub-haloes is only approximate. For exam- 
ple, host haloes of fixed mass show a fair amount of scatter 
in halo concentrations, correlated with the halo formation 
time (e.g., Wechsler et al. 2002; Zhao et al. 2003a;b; Lu 
et al. 2006), which we completely ignore. In addition, the 
retained mass fraction of sub-haloes is assumed to depend 
only on the instantaneous location of the sub-halo. In reality, 
however, sub-haloes with the same r a can have very different 
f m , depending on their orbital eccentricities and their time 
since being accreted by the host halo. Furthermore, eq. (10) 
was obtained for host haloes with masses ~ 10 14 /i _1 Mq, 
while we assume that it holds for haloes of all masses. Some 
of these shortcomings may affect our results. For example, 
as shown in Mandelbaum et al. (2005b), the scatter in the 
mass-luminosity relation can cause the derived mass to devi- 
ate from the mean sample halo mass, and so the uncertainty 
in the mass model of satellite galaxies may add uncertainty 
in our results. Unfortunately, a much more realistic model 
for the mass and density is not available at the present time, 
and we have to live with such uncertainties here. 



2.2 Central Galaxies 

We first consider the lensing signal around central galaxies. 
At sufficiently small R we may simply replace pCgm(»") in 
eqs (3) and (4) with the density distribution of the host 
halo, p(r). In what follows, we use AE c (i?) to refer to the 
ESD of central galaxies. 

For NFW profiles (i.e., a = l) we can write (see Wright 
& Brainerd 2000), 



ASc(i?) = 2^f J_1(C ' 1} [9{R/rc) ~ f{R/rc)] 



where 



(11) 



if x < 1 



/(*) = < 



and 




(12) 



Weak Lensing by Galaxies in Groups and Clusters 5 



100 - 



o 


w 50 



o 

< 



:>l i 1 i i i 1 i i i 1 i 
-■\\ 

i lrkoTKyf 1— 1 T n- 
-'. \\ !Og[M h J- 1 J.U, 


1 1 1 1 1 . 


- \ v 




- 


a = 0.0 


- i 


a=1.0 . 




a- 1 s 

LA — 1 . _ 




a = 2.0 _ 














_ ii 


i i 1 i i i 



0.2 0.4 0.6 
R (h^Mpc) 



0.8 



CM 
I 

o 

o 
2 



o 





-11 


i i i l i i i l i 




100 




!og[M h J— 1 J.O, 


OC— 1 .U _ 




i 




c = 3 








c = 6 




-i 




p- 1 


50 




c=15 





















i i i I 


* 1 . 


i i I i i i 







0.2 



0.4 
R (h- 



0.6 
»Mpc) 



0.8 



Figure 1. The excess surface density (ESD) around central galaxies as a function of radius R. Results are shown for 10 13 /i — 1 Mq dark 
matter haloes with different density profiles. The left-hand panel shows model predictions for haloes with a concentration parameter 
c = 10 but with different values for the central cusp slope, a. The right-hand panel shows the ESDs of NFW profiles (i.e., a = 1) with 
different concentration parameters, c. Note that in the presence of realistic measurement errors, there is a significant degeneracy between 
cusp slope and halo concentration. The arrows in both panels indicate the virial radius r v j r of a 10 13 ft -1 Mq dark matter halo. 
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g(x) = < 



As an illustration, Fig. 1 shows AE C as a function of 
the projected radius R for dark matter haloes with a mass 
Mh = 10 13 fe^ 1 M0. The left-hand panel depicts cases in 
which the halo concentration is fixed at c = 10, but in which 
the inner cusp slope a changes from to 2. The right-hand 
panel, on the other hand, shows results for NFW profiles 
(q = 1) but with different concentrations. As expected, for a 
given c, a model with smaller a has a shallower ESD at small 
radii. For the model with a constant density core (a = 0), 
the ESD goes to zero at the center, because there is no cen- 
tral density gradient that can cause image distortions. For a 
given a, more concentrated haloes have a more concentrated 
ESD, as expected (see also Guzik & Seljak 2002). Compar- 
ing the results shown in the left and right panels, one can see 
that there is a degeneracy between the central density gra- 
dient and the concentration of a dark matter halo, at least 
for realistic observational errors in the tangential shear. A 
similar degeneracy also hampers a unique derivation of the 
density distribution of dark matter haloes from the rotation 
curves of disk galaxies (e.g., van den Bosch et al. 2000). 



2.3 Satellite Galaxies 

The excess surface density around a satellite galaxy can be 
written as 



AX S (R\R S ) = AS s , sub (i?) + AS s , host (i?|i? s ) 



Here AE SjSU b(7?) and AE Sj host(i?|-Rs) are the ESDs due to 
the dark matter sub-halo and host halo, respectively, and 
R s is the projected distance between the satellite galaxy 

(13) 

and the center of its host halo (throughout we assume that 
a satellite resides at the center of its sub-halo). 

The azimuthally averaged, projected surface mass den- 
sity of the host halo around a satellite galaxy located at a 
projected distance R s from the halo center is 

Z s , bo st(R\R s ) = J E (y/R* + R? + 2R B RcosO) d0(15) 

where E(i?) is the projected density profile of the host halo. 
By integrating (15) from to R one obtains E Sj host(< R\R S ), 
and thus the ESD AE Sjho st(-R|^s)- The azimuthal averag- 
ing reflects the fact that in order to obtain sufficient signal- 
to-noise to measure the tangential shear one has to stack 
many satellite galaxies with the same R s - As long as these 
have random orientation angles 8, the azimuthal averaging 
of eq. (15) is appropriate. 

To illustrate the ESD around satellite galaxies, we 
first show the contribution from the host halo, by assum- 
ing that the satellite galaxy has no corresponding sub-halo 
(i.e., E SjSU b(-R) = 0). Fig. 2 shows the model predictions of 
AE Si host(^|^s) around satellite galaxies in NFW host haloes 
with masses of 10 12 h' 1 M Q (left panel) and 10 14 /i _1 M 
(right panel), respectively. Different line styles refer to dif- 
ferent halo-centric distances, as indicated. For comparison, 
we also show the results for R 3 = (solid lines), which 
are equivalent to the ESDs around central galaxies. Clearly, 
the ESD around satellite galaxies (without sub-haloes) are 
very different from those around central galaxies. They start 
with a value close to zero at R = 0, decrease to a neg- 
ative minimum near R s , and then increase rapidly with 
(14) radius, eventually approaching the ESD of central galax- 
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Figure 2. The excess surface density (ESD) around central galaxies (solid lines) and satellite galaxies as a function of radius R. Here it 
is assumed that satellite galaxies are not associated with dark matter sub-haloes; i.e., the ESD only reflects the contribution from the 
host halo. Different line styles correspond to different halo-centric distances of the satellite galaxies, as indicated. Panels of the left and 
right show the results for NFW host haloes with masses of 10 12 7i _1 M Q and 10 14 /i -1 M Q , respectively. 
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Figure 3. Same as Fig 2, except that this time satellite galaxies are assumed to reside in sub-haloes whose mass is equal to Af Sj o = 
10 ft -1 M0 at the time of accretion. The solid lines correspond to the ESD of a halo of this mass with R s = 0, and is shown for 
comparison. 



ies at R > 3R 3 ■ Note that here we are measuring the ESD 
AE s ,hostCR|-R s ) = E s , ho8t (< R\R S ) - E s>host (R\R S ) around 
satellite galaxies. When E Si host(-R|-Rs) reaches its maximum 
value at around R 3 , the ESD decreases to its negative mini- 
mum. As expected, the overall amplitude of AE 3j host(ii|iis) 
at large radii is higher for more massive haloes. Thus, if the 
value of R s is known, the tangential shear 7t(i?) measured 
around satellite galaxies can be used to constrain the mass 



distribution of their host haloes. We will come back to this 
in Section 3.2. 

We now include the contribution of the sub-halo. For 
this, we use the sub-halo model described in Section 2.1, 
which allows us to compute AE SjSU b(-R) for a given sub- 
halo mass. This mass can be related to the sub-halo mass 
at the time of accretion, using eq. (10) and the true dis- 
tance r s of the sub-halo from the center of the host halo. 
In order to compute the contribution due to the host halo, 
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AS Sj host(^|^s), one also needs to know R s , which is the 
projection of r s on the plane of the sky (i.e., < R s < r s ). 

As an illustration, we consider a simple case in which 
all sub-haloes have an original mass M Sj0 = 10 ll s /i -1 M Q 
at the time of accretion. We model the tidal mass loss as 
described in Section 2.1 and adjust their density profiles ac- 
cordingly. Fig 3 shows the full ESD (including both the host 
and the sub-halo terms) around satellite galaxies located 
at different halo-centric distances in host haloes of different 
masses. For simplicity we assume that R s — r s , i.e., that the 
radius vector f s is perpendicular to the line-of-sight. Com- 
paring Figs 3 and 2, one sees that the contribution of the 
sub-halo to the ESD dominates at small R, as expected. As 
a comparison, the solid curves in Fig 3 show the AE SjSU b(i?) 
of sub-haloes with a mass M s = 10 ll s /i _1 Mq (the contri- 
bution from the host-term is ignored here). At small radii 
these are slightly higher than those represented by the dot- 
ted and dashed curves, because sub-haloes are modeled to 
have lost a fraction 1 — f m of their mass after having been 
accreted by the host halo. Note that satellites with smaller 
R s have a somewhat lower ESD at small radii. This owes 
to the fact that the mass loss is larger for sub-haloes at a 
smaller halo-centric distance. 



3 TEST USING GALAXIES AND GROUPS IN 
MOCK CATALOGUES 

The analysis presented above shows that it is in principle 
possible to use galaxy-galaxy lensing to probe the mass dis- 
tribution of haloes and sub-haloes. Unfortunately, we do not 
know a priori whether a galaxy is a central galaxy or a satel- 
lite galaxy. Since the expected galaxy-galaxy lensing signal 
of individual galaxies is very weak, one has to combine the 
shear measured around a large number of galaxies to ob- 
tain a statistically significant detection. The ESD inferred 
from such a measurement is therefore the average over all 
lensing galaxies used. It is clear from the above, that if that 
average combines both central and satellite galaxies, the re- 
sulting AE(_R) will be difficult to interpret, especially when 
the satellites have different R s . 

However, if we could identify central and satellite galax- 
ies a priori, and if we could also determine R s for each satel- 
lite, we could measure the tangential shear of central and 
satellite galaxies separately, thus constraining the mean host 
halo mass as well as the mean sub-halo mass of the satel- 
lites directly. Such an identification of central and satellite 
galaxies requires a method to decide which galaxies belong 
to the same dark matter halo. In a recent paper, Yang et 
al. (2005a) developed a new, halo-based group finder which 
is particularly successful in grouping galaxies in a redshift 
survey according to their common dark matter haloes. In 
this section, we apply this group finder to mock redshift cat- 
alogues and examine whether the membership information 
thus obtained is useful in galaxy-galaxy lensing studies. In 
particular, we examine whether the identification of central 
and satellite galaxies inferred from such a group catalogue is 
sufficiently reliable that it allows an accurate measurement 
of the masses of the corresponding host- and sub-haloes. 



3.1 N-body simulations and mock catalogues 

We use the results of a high-resolution N-body simula- 
tion to construct mock galaxy catalogues. The simulation 
was carried out on the VPP5000 Fujitsu supercomputer 
of the National Astronomical Observatory of Japan with 
the vectorized-parallel P 3 M code (Jing & Suto 2002). It 
evolves the distribution of the dark matter from an ini- 
tial redshift of z = 72 down to z — in a ACDM 'con- 
cordance' cosmology with £l m — 0.3, Qa = 0.7, h = 
-ffo/(100 kms -1 Mpc -1 ) = 0.7 and with a scale-invariant 
initial power spectrum with normalization as — 0.9. The 
simulation uses 512 3 cold dark matter particles in a peri- 
odic cube of 100 x 100 x W0h~ 3 Mpc 3 . The particle mass is 
equal to 6.2 x lO 8 /i -1 M0, and dark matter haloes are identi- 
fied using the standard FOF algorithm with a linking length 
of 0.2 times the mean inter-particle separation. 

In order to construct realistic mock galaxy samples, 
we populate the haloes with galaxies, using the conditional 
luminosity function (hereafter CLF; Yang, Mo & van den 
Bosch 2003b; van den Bosch, Yang & Mo 2003), $(L|M), 
which gives the average number of galaxies of luminosity L 
that reside in a halo of mass M. The CLF model parame- 
ters used to construct the mock galaxy catalogue are given 
in Table 1 of van den Bosch et al. (2005c) as ID # 6. We re- 
fer the reader to Yang et al. (2004) for details regarding the 
construction of these mock galaxy catalogues. We emphasize 
though, that by construction, these mock catalogues match 
the observed luminosity function and the observed cluster- 
ing strength as function of luminosity. 

The CLF model allows us to populate each halo in the 
simulation with galaxies of different luminosities. We locate 
the brightest galaxy in each halo at the halo center and 
assume that the other galaxies (the satellites) are associ- 
ated with dark matter sub-haloes. Unfortunately, we can 
not use the actual sub-haloes in the numerical simulation 
itself, simply because the resolution of the simulation is not 
sufficient to resolve sub-haloes in low mass host haloes. In 
addition, the survival and structure of sub-haloes is affected 
by the baryonic component, which is not modeled in our 
dark-matter-only simulation. We therefore follow a different 
approach. We distribute the satellite galaxies isotropically 
throughout the host halo, with a number density distribu- 
tion that reproduces the dark matter density profile. Note 
that here we use the spherical NFW profile with the con- 
centration obtained by Eke, Navarro & Steinmetz (2001) in 
modeling the spatial distribution of the satellite galaxies rel- 
ative to the halo center. Next, we model the sub-halo popu- 
lation associated with these satellite galaxies adopting a pre- 
scription similar to that used in Kravtsov et al. (2004), Vale 
& Ostriker (2006), and Conroy et al. (2006): we assume that 
the luminosity of the satellite has a monotonic relation with 
the sub-halo mass at the time of its accretion. In order to es- 
tablish such a relation, we use the mass function of sub-halo 
progenitors obtained by van den Bosch et al. (2005a), shown 
in the left-hand panel of Fig. 4, and the conditional luminos- 
ity function. The sub-halo mass function has been studied 
recently both in numerical simulations and galaxy-galaxy 
lensing observations (e.g. De Lucia et al. 2004; Natarajan & 
Springel 2004; Natarajan et al. 2006). The mass- luminosity 
relations thus obtained are shown in the right-hand panel 
of Fig. 4, for four different host halo masses. Using these 
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Figure 4. The left-hand panel shows the un-evolved mass function of dark matter sub-haloes from van den Bosch et al. (2005a). Here 
'un-evolved' means that the masses correspond to the masses of the sub-haloes at the time that they were accreted (i.e., before any mass 
loss due to tidal stripping occurred). As shown in van den Bosch et al. (2005a), this mass function is independent of the mass of the host 
halo. The right-hand panel shows the relation between the luminosity of a satellite galaxy and the mass of the sub-halo at the moment 
of accretion. This relation is obtained by comparing the number density of satellite galaxies predicted by the conditional luminosity 
function with the number density of sub-haloes predicted by the un-evolved sub-halo mass function shown in the left-hand panel. 



luminosity-mass relations and the model of sub-halo struc- 
ture described in Section 2.1, we assign each satellite a sub- 
halo mass depending on (i) its luminosity, (ii) its host halo 
mass, and (iii) its distance from the center of the host halo. 

Next we construct a mock redshift catalogue. We put 
the center of the simulation box at a comoving distance of 
200 /i _1 Mpc from a virtual observer. Each galaxy is given 
a redshift according to its distance from this observer and 
its peculiar velocity along the corresponding line of sight. 
We construct a flux-limited sample by including only those 
galaxies that have an apparent magnitude nib, < 19.30. Fi- 
nally, we apply our halo-based group finder to this mock 
galaxy catalogue in order to construct a catalogue of galaxy 
groups. Following Yang et al. (2005b), we assign a halo mass 
to each group according to its ranking in total group lumi- 
nosity and using the halo mass function for the standard 
ACDM cosmology (see Yang et al. 2005b for more details). 
Tests in Yang et al. (2005b) and Weinmann et al. (2006) have 
shown that, on average, the group masses thus obtained are 
in good agreement with the input halo masses. 

In what follows we refer to the brightest galaxy in each 
group as the central galaxy, and to all other group members 
as satellites. We use the mock galaxy group catalogue to 
test how well galaxy-galaxy lensing measurements around 
these 'central' and 'satellite' galaxies allow a recovery of 
their host and sub-halo masses respectively. Due to inter- 
lopers (group members that do not belong to the same dark 
matter halo) and the fact that the group finder may occa- 
sionally miss a halo member, central and satellite galaxies 
in the group catalogue are not necessarily also central and 
satellite galaxies in their real haloes. If this confusion is too 
large, it will not be possible to obtain reliable estimates of 
host and/or sub-halo masses. Clearly, the accuracy of such 



an approach therefore needs to be tested, which is the pur- 
pose of the mocks constructed here. In the left-hand panel 
of Fig. 5, we show the fractions of the false central galax- 
ies and false satellite galaxies (interlopers) as a function of 
group mass. The fraction of false central galaxies is com- 
pletely negligible in groups of all masses, while that of false 
satellites is typically well below 20%. As we have pointed out 
in section 2.3, the ESDs around satellite galaxies depend 
strongly on the projected halo-centric distances. Therefore 
it is also important to check how the interloper fractions of 
the satellite galaxies depend on the projected group-centric 
distances. The results are shown in the right-hand panel of 
Fig. 5. One may notice that at small projected group-centric 
distances the interloper fractions are rather small, and then 
increase (due to the larger area covered) as the projected 
group-centric distances increase. We will discuss later how 
and to what extent this interloper fraction may affect our 
measurements of the ESDs and the extracted properties of 
the subhaloes. 



3.2 Lensing by galaxies in groups 

3.2.1 The lensing signal around central galaxies 

We first consider the lensing signal around the central group 
galaxies. In order to "measure" the ESD we project the po- 
sitions of galaxies and dark matter particles onto a plane 
perpendicular to the line of sight, and estimate the mean 
dark matter surface density contrast within rings of differ- 
ent radii around the galaxies. Since by definition the back- 
ground surface density is subtracted, the ESD thus derived 
is independent of the depth of the projection (here equal to 
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Figure 5. The left-hand panel shows the fractions of the false central (solid line) and satellite galaxies (dashed line) as a function of 
group masses. The right-hand panel shows the fractions of the false satellite galaxies (interlopers) in groups of different mass bins (as 
indicated) as a function of the projected group-centric distances. 
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Figure 6. The ESDs measured from the N-body simulation around central galaxies (open circles) and around satellite galaxies at different 
R s , as indicated (open triangles and solid squares) Note that here satellite galaxies are modeled without dark matter sub- haloes. The 
two panels correspond to different bins of assigned group mass, as indicated with masses in h _1 Mq. The lines in each panel are the 
best-fit models to the ESDs shown: For central galaxies, this model has two free parameter (M^ and c), while in the case of satellite 
galaxies it has four free parameters (M^, c, M s and R s )- Since we do not have realistic 'measurement' errors, we give each 'data point' 
equal weight in the fitting. See text for more details. 



the simulation box size), as long as this depth is much larger 
than the dark matter correlation length. 

The open circles in Fig 6 show the AE(i?) thus obtained 
around the central galaxies in groups with assigned masses 
in the range 12.0 < log(M h //i _1 M©) < 12.5 (left panel) and 
14.0 < log(Afh/ 7i _1 M ) < 14.5 (right panel), respectively. 
These ESDs are monotonic functions of radius, and the am- 



plitude is higher for the more massive haloes, in agreement 
with the predictions shown in Section 2. In Section 3.3, we 
investigate the accuracy with which these measurements al- 
low a recovery of the mean host halo mass of these galaxies. 
Note that the group masses have been assigned based on 
the total group luminosity; the lensing signal will provide a 
direct test of this mass. 
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3.2.2 The lensing signal around satellite galaxies 

As discussed in Section 2, the galaxy-galaxy lensing signal 
from satellite galaxies depends not only on the luminosity of 
the satellite, but also on the properties of the host halo and 
the halo-centric distance. This means that the tangential 
shear measured around a large sample of satellite galaxies in 
a given luminosity bin is very difficult to interpret. Rather, 
one can use the actual information from the group catalogue 
to sort the satellites according to both the mass of their 
host halo (i.e., the assigned group mass) and the projected 
distance from the center of the group. Such a signal is much 
easier to interpret, as it tightly constrains the contribution 
of the AE Sj host-term to the shear measurements. 

To illustrate the potential power of this approach we 
measure the ESD around satellite galaxies in our mock group 
catalog, using the projected distribution of dark matter par- 
ticles as in Section 3.2.1 above. Since the simulations do not 
resolve the majority of sub-haloes, and since we did not as- 
sociate satellite galaxies in the mock with the resolved sub- 
haloes, this ESD reflects the lensing signal that one would 
obtain if satellite galaxies are not surrounded by dark matter 
sub-haloes. 

The open triangles and solid squares in Fig 6 show 
the AS(i?) around satellite galaxies with different projected 
distances from their central group galaxies, as indicated. 
As demonstrated in Section 2.3, the ESDs around satel- 
lite galaxies are expected to become similar to those around 
their central galaxies at radii R>3R S . However, in Fig 6 
the ESDs of the satellite galaxies have slightly higher ampli- 
tudes at large radii than those of the corresponding central 
galaxies. This systematic offset is due to the fact that wc 
have combined groups in a finite mass range. Since more 
massive haloes (groups) in general contain a larger number 
of satellite galaxies, any satellite-averaged mean, such as the 
ESD, will be biased towards the more massive haloes. A sim- 
ilar bias occurs when one tries to estimate the average halo 
mass of a stack of host galaxies from the velocity dispersion 
of their satellite galaxies (see van den Bosch et al. 2004) . 

The ESD around satellite galaxies in groups with 
10 12 /i _1 M Q < M h < 10 12 ' 5 /i _1 M Q and with R s = 
0.2 ± 0.05k- 1 Mpc (solid squares in left-hand panel) reveals 
a small upturn at i?<0.1^ _1 Mpc. However, according to 
Fig. 2, AE(J?) should go to zero at small R. This disagree- 
ment with the theoretical predictions is caused by the fact 
that a small number of these satellite galaxies in the group 
catalogue are interlopers, which in reality are mainly central 
galaxies in lower mass haloes. Since the interloper fraction 
increases with halo-centric distance (shown in the right-hand 
panel of Fig. 5), this explains why this peak near R = is 
more pronounced for satellites with a larger R 3 . Clearly, it is 
important that any group finder used for identifying centrals 
and satellites is properly calibrated to yield sufficiently low 
interloper fractions. As described in Yang et al. (2005a), the 
group finder used here has been calibrated accordingly. In 
particular, the interloper fraction is constant (at ~ 20 per- 
cent) with group mass, unlike the more standard friends-of- 
friends method which typically yields an interloper fraction 
that increases systematically with decreasing group mass 
(see Fig. 7 in Yang et al. 2005a). 

We now add the contribution of the dark matter sub- 
haloes. The density distribution of a sub-halo surrounding a 



satellite galaxy of a given luminosity, located at a given halo- 
centric distance in a host halo of a given mass is modeled as 
described in Section 3.1, and its contribution to the lensing 
signal is computed by integrating p s (r) along the line-of- 
sight. The resulting ESDs are shown in Fig. 7. As we have 
seen before, the contribution of the dark matter sub-haloes 
enhances the ESD around satellite galaxies on small scales. 

3.3 Estimating the masses of haloes and 
sub-haloes from galaxy-galaxy lensing 

Having estimated the lensing signal that one could in princi- 
ple measure around central and satellite galaxies (given vir- 
tually infinite signal-to- noise in the shear measurements), 
we now investigate to what extent such a signal allows to 
recover the average masses and concentration parameters of 
both the host haloes and the sub-haloes. 



3. 3. 1 The masses of haloes 

We start by testing how accurately the lensing signal around 
galaxies classified as 'centrals' by our group-finder allows a 
recovery of the mean mass and concentration of their host 
haloes. To that extent we use the assigned group masses as 
a 'pre-selection', and compute the AS c (i?) around central 
galaxies in three group- mass bins, two of which are shown in 
Fig 6. Assuming that the density distributions of dark mat- 
ter haloes are well described by NFW profiles, we fit these 
AS c (i?) using eq. (11) with c and Mh as free parameters. 
The best fits are shown as solid lines in Fig. 6, and should 
be compared to the open circles. 

The corresponding best fit values of Mh and c are shown 
in the left panel of Fig 8 as the open circles. For comparison, 
the asterisks show the true mean halo masses that host the 
mock galaxies used to determine AE C (_R), while the hatched 
area indicates the bin of assigned group masses. Note that 
the best-fit Mh is in good agreement with the true average, 
and that both lie well within the range of assigned group 
masses. The recovered concentration parameters also match 
the expected values predicted by Eke, Navarro & Stein- 
metz (2001), indicated by the solid line, reasonably well. 
Although there are slight deviations, they are significantly 
smaller than the typical scatter in c for haloes of a given 
mass, Alogc ~ 0.14, indicated by the errorbars attached to 
the asterisks (e.g. Jing 2000, Bullock et al. 2001, Wechsler 
et al. 2002). 

In the case of real galaxy-galaxy lensing data, the errors 
on AS c (i?) will depend on the details of the selection of the 
lens-source samples and on the quality of the observational 
data. We do not model such errors here; instead we give 
each data point an equal weight in the fitting. We have also 
tested the impact of changing the weight scheme, the change 
in the best fitting parameters are very small. Note that the 
errorbars from recent galaxy-galaxy lensing measurements 
show different radial dependences (e.g. Mckay et al. 2001; 
Mandelbaum et al. 2005a), the absence of realistic errorbars 
on our 'measurements' prevents us from putting meaningful 
confidence levels on the best-fit parameters. 

The ESDs measured around satellite galaxies also con- 
tain information regarding the mass of the host halo. In 
order to see how well these lensing measurements allow to 
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Figure 7. Same as Fig. 6, except that here satellite galaxies are associated with dark matter sub-haloes, as described in the text. Note 
that this enhances the ESD at small R. 



recover the masses of their host and sub-haloes we fit with 
eq. (14) the AE S (R\R 3 ) for a number of bins in R s . We have 
four free parameters in the fit: Mh, c, M s and R s (subscripts 
h and s refer to the host and sub-haloes, respectively). Note 
that c is the concentration of the host halo, not that of 
the sub-halo. Since the inner part of the ESDs cannot be 
measured with high precision, we assume that all sub-haloes 
have a concentration c 3 — 10. We will discuss later in section 
3.3.2 how this assumption may affect the extracted proper- 
ties (i.e. masses) of the sub-haloes. Note also that we let R a 
be a free parameter, eventhough we only consider satellite 
galaxies in a relatively narrow bin in R„ (as determined from 
the group catalogue). 

The best-fit curves are given by the solid lines in Figs. 6 
and 7. Note that these best-fit models are in reasonable 
agreement with the 'measured' ESD. The slight deviations 
are due to interlopers (i.e., errors associated with the group 
catalogue) and to the fact that the model does not account 
for the fact that the data corresponds to a range in host 
and sub-halo masses (i.e., AE S (_R|_R S ) reflects a linear com- 
bination of many slightly different ESDs). The best-fit host 
halo masses are shown in the right panel of Fig. 8: solid 
dots are the results obtained using the mock in which we 
do not include the analytical sub-haloes (Fig. 6), while the 
open circles show the results obtained when sub-haloes are 
included (Fig. 7). The agreement with the true mean host 
halo masses, indicated by the asterisks, is extremely good, 
suggesting that the lensing signal from satellite galaxies ba- 
sically yields an equally good measurement of the host halo 
mass as that from the central galaxies. Thus, with a well- 
defined group catalog, the properties of the host haloes can 
be accurately recovered from the ESDs using both central 
and satellite galaxies. 

There is one small caveat though. Note that the host 
halo masses obtained from the central galaxies are some- 
what lower than those obtained from the satellites. This is 
again due to the fact that the ESDs of satellite galaxies are 



weighted by the number of satellite galaxies; since there are 
more satellite galaxies in more massive groups, more massive 
systems receive a larger weight in the averaging. In partic- 
ular, groups that consist only of a single member (which 
then by definition is a central galaxy), only contribute to 
AE C , but not to AE S . Since such a group is likely to have 
a relatively low mass compared to the other groups in the 
same mass bin, this will bias the value of Mh inferred from 
the satellites high with respect to the value of Mh inferred 
from the central galaxies. Note, however, that this does not 
reflect an error; the average true halo mass of systems with 
at least one satellite is simply larger than the average true 
halo mass of all systems (including those with zero satel- 
lites). This is also evident from the fact that the asterisks 
in the left and right hand panels of Fig. 8 indicate differ- 
ent masses. Thus, if the host masses inferred from AE C and 
AE S do not agree with each other, this does not necessar- 
ily indicate an inconsistency. It may also simply reflect a 
'selection effect', in that the ESDs around central and satel- 
lite galaxies are contributed by somewhat different host halo 
populations. As a simple test, we weight the ESDs for cen- 
tral galaxies with the number of satellite galaxies at different 
group-centric distances. The recovered host halo masses are 
now in much better consistent with those recovered from the 
ESDs around satellite galaxies. 



3.3.2 The masses of sub-haloes 

The four-parameter fits to the ESDs of satellite galaxies also 
yield best-fit values for the average sub-halo mass, M s and 
for the projected distance, R s . The open circles in Fig. 9 
show the best-fit values of M s and R s as obtained from the 
mock in which we added analytical sub-haloes around each 
satellite galaxy (i.e., these correspond to the ESDs shown 
in Fig. 7). For comparison, the asterisks show the corre- 
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Figure 8. Properties of host haloes recovered from fitting the ESDs around central galaxies (left-hand panel) and satellite galaxies 
(right-hand panel) in the mock group catalogue. The hatched areas in both panels indicate the range of assigned masses of the groups 
from which the central and satellite galaxies are selected. The open circles in the left-hand panel indicate the masses and concentrations 
of the NFW haloes that best-fit the ESDs of the central group galaxies. The asterisks show the true mean host halo masses and the 
corresponding halo concentrations obtained using the model of Eke et al. (2001; shown as solid line). The errorbars on the asterisks 
correspond to A log c = 0.14 which reflects the typical scatter in c at a given halo mass. The solid dots and open circles in the right-hand 
panel indicate the best-fit host halo masses and group-centric radii obtained from fitting the ESDs around satellite galaxies with and 
without sub-halo contributions, respectively. The horizontal error bars on these points indicate the ranges of group-centric radius, R a , 
that were used to select the satellite galaxies from the group catalogue. The asterisks indicate the true mean host halo masses and the 
true mean R a , for the corresponding satellite galaxies. 
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Figure 9. Open circles indicate the best-fit sub-halo masses and group-centric radii obtained from fitting the ESDs around satellite 
galaxies in the mock group catalogue in which each satellite galaxy has been assigned a dark matter sub-halo as described in Section 3. 
Different panels correspond to groups in different bins of assigned mass, as indicated, while the horizontal error bars indicate the ranges 
of R B that were used to select the satellite galaxies from the group catalogue. The hatched areas show the 90 percentile ranges of true 
(retained) sub-halo masses. The asterisks indicate the true mean sub-halo masses and the true mean R s , for the corresponding satellite 
galaxies. Finally, the solid squares (also with horizontal errorbars) indicate the best-fit values of M s and R a obtained by fitting the ESDs 
around the satellite galaxies in the mock without dark matter sub-haloes. The fact that these best-fit values are not equal to zero, as 
they should be, owes to the contamination due to interlopers. See text for a detailed discussion. 
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Figure 10. Same as Fig. 9, except that here we show the results for different bins in satellite luminosity L s , as indicated in h~ 2 Lq. 
As in Fig. 9 the open circles show the best-fit sub-halo masses and group-centric radii obtained fitting the ESDs assuming NFW sub- 
haloes with c s = fO. The solid diamonds show the same quantities but obtained by assuming sub-haloes modeled by eq. (7) with mean 
<r s ) « (Rs +r vir , h )/2. 



sponding true mean values 1 . The solid squares in Fig. 9 
indicate the best-fit values of M 3 obtained from the mock 
in which no sub-haloes are included (i.e., corresponding to 
the ESDs shown in Fig. 6). The fact that these best-fit val- 
ues are not equal to zero, as they should be, owes to the 
presence of interlopers that in reality are centrals. These 
best-fit sub-halo masses should therefore be considered a 
contamination of the measurement in the case with sub- 
haloes. Typically this contamination is much smaller than 
the actual mean sub-halo mass, especially at small R 3 . In 
the cases examined here, the largest contamination we find 
is for the satellites with R s = 0.3 ± 0.05/i~ Mpc in groups 
with 10 12 /i _1 M Q < M h < 10 12 ' 5 /i _1 M s , where the best-fit 
M a inferred from the mock without sub-haloes is 38 per- 
cent of that obtained with sub-haloes. In all other cases this 
fraction is much lower. 

The shaded areas in Fig. 9 outline the 90 percentile 
range of the real sub-halo masses, which clearly is very 
broad. This owes to the fact that we have added the lensing 
signal from all satellite galaxies, irrespective of their lumi- 
nosity. In our model, and most likely also in reality, more 
luminous satellites reside in more massive sub-haloes. This 

1 For interlopers that are centrals, rather than satellites, we use 
the host-halo mass in the computation of the averages. 



suggests that the actual range of sub-haloes probes may be 
narrowed by selecting satellite galaxies of similar luminosi- 
ties. Fig. 10 shows the ESDs around satellite galaxies in dif- 
ferent luminosity bins and in groups of different (assigned) 
masses. As in Fig. 9, the shaded areas indicate the 90 per- 
centile ranges of the true sub-halo masses. As expected, 
these are now much narrower. The remaining width owes 
to the finite widths of the ranges in satellite luminosity, L s , 
group mass, Mh, and projected separation, R s , and to the 
fact that projection causes systems with different r s to con- 
tribute to the same R s ■ Note that in our model the sub-halo 
mass is a function of Mh, L s and r s . 

The open circles in Fig. 10 indicate the best-fit sub- 
halo masses, obtained by fitting the corresponding ESDs 
with a four parameter model as described above. Note that 
the recovered sub-halo masses are slightly larger than the 
input value, especially for fainter satellite galaxies in massive 
groups. We propose that the main reason for this is that a 
stripped halo has a steeper mass profile (see eq. 7), and 
the assumption of a fixed NFW profile with c 3 — 10 is not 
sufficiently accurate. As a simple test of this hypothesis, we 
consider a simple model in which we assume that the mean 
halo-centric distances of the satellite galaxies in each sample 
is (r s ) « (R s + r v i T> h)/2. Then, we use eq. 7 to model the 
mean profile of the sub-haloes in consideration. This model 
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ensures that sub-haloes with smaller R 3 , which are expected 
to have suffered more mass loss, are more concentrated. The 
best-fit sub-halo masses obtained from this model are shown 
in Fig. 10 as the solid diamonds. As can be seen, the best-fit 
masses are systematically smaller than in the case of c 3 = 10, 
and are in better agreement with the true averages. 

In conclusion, with a well-defined group catalog, the 
ESDs of satellite galaxies obtained from galaxy-galaxy lens- 
ing, with signals properly stacked according to the group 
masses (inferred from the total group luminosities) and the 
luminosities and group-centric distances of satellite galaxies, 
can be used to probe the masses of sub-haloes associated 
with satellite galaxies. 

3.3.3 What if the brightest galaxy is off-center? 

In the analysis presented above, the brightest halo galaxy is 
assumed to reside at the halo center. In reality, the location 
of the brightest galaxy may not be exactly at the halo center, 
and this may affect the predicted ESDs around the central 
and satellite galaxies. Here we examine the importance of 
this effect. 

In a recent study of the phase space distribution of 
the brightest halo galaxies in galaxy groups, van den Bosch 
ct al. (2005b) found that, for haloes with masses Mh > 
10 13 /i _1 Mq, the deviation of the location of the brightest 
galaxy from the halo center is about 3 percent of the halo 
virial radius. Based on SPH simulations, Berlind (2003) also 
found that the 'central' galaxy in a dark matter halo may 
deviate from the position of the most bound particle by 2 
percent of the halo virial radius. In order to quantify how 
such deviations impact on the ESDs, we perform the follow- 
ing test. We construct a new mock catalogue in which we 
assume that the position of each brightest halo galaxy devi- 
ates from the most bound particle with an amount given by 
a Gaussian distribution with a dispersion equal to 3 percent 
of the halo virial radius. Next we run our group finder over 
this new mock catalogue, compute the ESDs around the cen- 
tral and satellites galaxies, and fit them as described above 
to determine the best-fit halo masses and concentrations. 
We find that the off-centering of the central galaxies has an 
almost negligible effect, except that the best-fit concentra- 
tion parameters of host haloes are reduced by a modest 15 
percent. Such a reduction arises from the fact that the off- 
centering decreases the ESD at small radius, similar to what 
happens to the ESD around satellite galaxies. Thus, the abil- 
ity to measure the host halo masses is not compromised by 
the fact that central galaxies may be (slightly) off-centered 
from the gravitational center of their dark matter halo. 



4 DISCUSSION AND CONCLUSIONS 

The weak galaxy-galaxy lensing signals from central and 
satellite galaxies (at different halo-centric distances) are very 
different. So far, this difference has not been fully exploited. 
In fact, in all previous studies it has been considered a nui- 
sance, rather than a source of valuable information. Most 
galaxy-galaxy lensing analyses have either simply ignored 
the differences between central and satellite galaxies, or have 
focused on isolated galaxies in an attempt to minimize the 
disturbing contribution of satellite galaxies. An alternative 



approach, pioneered by Guzik & Seljak (2002), is to model 
the combined lensing signal from central and satellites. How- 
ever, no study to date has attempted to measure the lensing 
signals of centrals and satellites at different halo-centric dis- 
tances separately. This simply owes to the fact that it is not 
a priori clear which galaxy is a central galaxy and which 
galaxy a satellite. This is unfortunate, as a separation of 
these two components would allow a much cleaner measure- 
ment of the halo masses (and their detailed density profiles) 
hosting central galaxies. In addition, it would allow for a 
direct measurement of the masses (and density profiles) of 
sub-haloes hosting satellite galaxies. Furthermore, since the 
shear around satellite galaxies also harbors information re- 
garding the host haloes in which their sub-haloes reside, one 
can use their lensing signal to obtain an independent mea- 
surement of the host halo masses. 

In order to exploit this richness in information, it is cru- 
cial that one has a reliable technique for separating central 
and satellite galaxies. In this paper, we have investigated 
to what extent this is feasible with the halo-based group 
finder developed by Yang et al. (2005a) . This group finder is 
straightforward to apply to galaxy redshift surveys, and has 
been calibrated to yield low interloper fractions. The lat- 
ter is important to minimize the confusion of centrals and 
satellites. In order to test our methodology, we constructed 
detailed mock redshift surveys from N-body simulations that 
are populated with galaxies using the conditional luminos- 
ity function. Application of the halo-based group finder to 
this mock redshift survey yields a large catalogue of mock 
galaxy groups, from which central galaxies are identified as 
the brightest galaxy in each group. All other group members 
are considered to be satellites. 

The group catalogue is not only useful to separate cen- 
trals from satellites. It also yields, for each satellite, the 
projected separation, R s , to it central galaxy. This is ex- 
tremely important, since the contribution of the host halo 
to the gravitational shear around a satellite galaxy depends 
strongly on R s . Therefore, in order to facilitate a meaning- 
ful interpretation of the shear measurements around satellite 
galaxies, it is crucial that one only stacks the data from a 
relatively narrow bin in R s . Another piece of useful infor- 
mation from the group finder are the masses that one can 
assign to each group based on the total luminosity of all 
its members. This mass can be used to 'pre-select' central 
and satellite galaxies in relatively narrow bins in group mass, 
which again facilitates the interpretation of the lensing mea- 
surements. 

To test these ideas we computed the excess surface den- 
sity (ESD), which is the observable that can be obtained 
from the lensing measurements, around each central and 
satellite galaxy in the ./V-body simulation. We have demon- 
strated that by stacking the ESDs around central galaxies 
within different group mass bins, the average masses and 
concentrations of their host haloes can be accurately recov- 
ered from the data. In addition, we have shown that the 
ESDs around satellite galaxies, when stacked according to 
group-centric distance and group mass, allow an equally ac- 
curate recovery of the masses of their corresponding sub- 
haloes, as well as that of their host haloes. 

We therefore conclude that a combination of galaxy- 
galaxy lensing measurements with a galaxy group cata- 
logue extracted from a large redshift survey, such as the 
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Sloan Digital Sky Survey, in principle allows for accurate 
measurements of the masses and concentrations of host 
haloes around central galaxies and sub-haloes around satel- 
lite galaxies. However, it is important to realize that we have 
not attempted to mimic realistic observations. Rather, we 
have simply assumed infinite accuracy in the measurements 
of the ESDs. In reality, resolution issues due to the finite 
sampling of the shear field and errors in the shear measure- 
ments, may cause a substantial reduction of the accuracy 
with which this methodology can be applied. In a forth- 
coming paper, we will apply this method to realistic mock 
shear maps taking into account various observational effects 
to test its feasibility. 
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